Dynamics of fluctuations in a fluid below the onset of Rayleigh-Benard convection 
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We present experimental data and their theoretical interpretation for the decay rates of temper- 
ature fluctuations in a thin layer of a fluid heated from below and confined between parallel hori- 
zontal plates. The measurements were made with the mean temperature of the layer corresponding 
to the critical isochore of sulfur hexafluoride above but near the critical point where fluctuations 
are exceptionally strong. They cover a wide range of temperature gradients below the onset of 
Rayleigh-Benard convection, and span wave numbers on both sides of the critical value for this on- 
set. The decay rates were determined from experimental shadowgraph images of the fluctuations at 
several camera exposure times. We present a theoretical expression for an exposure-time-dependent 
structure factor which is needed for the data analysis. As the onset of convection is approached, 
the data reveal the critical slowing-down associated with the bifurcation. Theoretical predictions 
for the decay rates as a function of the wave number and temperature gradient are presented and 
compared with the experimental data. Quantitative agreement is obtained if allowance is made 
for some uncertainty in the small spacing between the plates, and when an empirical estimate is 
employed for the influence of symmetric deviations from the Oberbeck-Boussinesq approximation 
which are to be expected in a fluid with its density at the mean temperature located on the critical 
isochore. 
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I. INTRODUCTION 

In this paper we present experimental data for and 
a theoretical analysis of the decay rates of fluctuations 
in a fluid layer between two horizontal plates that are 
maintained at two different temperatures. The size L 
of the layer in the horizontal directions is much larger 
than the distance d between the plates. By now it has 
been well established that the presence of a uniform and 
stationary temperature gradient VTq in the fluid layer 
induces hydrodynamic fluctuations that are long ranged 
in space. 

An important dimensionless parameter that governs 
the nature of the thermal noncquilibrium fluctuations is 
the Rayleigh number 



R = 



agd-^VTo 



(1) 



where g is the gravitational acceleration and where a is 
the isobaric thermal expansion coefflcient, Dt the ther- 
mal diffusivity, and v the kinematic viscosity of the fluid. 
The Rayleigh number is commonly taken to be positive 
when the fluid layer is heated from below. States in which 
the fluid layer is heated from above then correspond to 
i? < 0. The fluid layer in the presence of a tempera- 
ture gradient remains in a quiescent state for all R less 
than a critical value R^ for the onset of Rayleigh-Benard 
(RB) convection, including all negative values of R. The 
present paper is concerned with thermal fluctuations in 
such a fluid layer for positive values of R below Re- 

At sufflciently large wave numbers q the inten- 
sity of noncquilibrium fluctuations is proportional to 
(Vro)2 g-4 both for negative and positive i? 0,111111. 



For smaller values of q it was shown theoretically |1| and 
confirmed experimentally 6] that the increase of the fluc- 
tuation intensity with decreasing q is quenched in the 
presence of gravity, yielding a constant value in the limit 
of small q. If the presence of top and bottom boundaries 
is taken into account, one finds that the intensity of the 
fluctuations vanishes as q^ '3', '3| at small q. Hence, 
the nonequilibrium structure factor S{q) is predicted to 
exhibit a crossover from a dependence for larger q to 
a q^ dependence in the limit g — > 0, leading to a max- 
imum at an intermediate wave number (/max which has 
a value near iv/d. As the RB instability is approached 
from below, linear theory predicts that S'((7max), as well 
as the total fluctuation power [i.e., the integral of S{q)\ 
diverges, in agreement with the asymptotic prediction 
obtained by Zaitsev and Shliomis [9j and by Swift and 
Hohenberg [HllJ. The predicted increase as the RB in- 
stability is approached was verified quantitatively by ex- 
periments using shadowgraph measurements . As the 
fluctuations become large it was predicted |lOl Hlf and 
confirmed experimentally that linear theory breaks 
down and that the fluctuation amplitudes saturate due 
to nonlinear interactions. 

The present paper is concerned not with the inten- 
sity, but with the dynamics of the nonequilibrium fluc- 
tuations. One experimental technique to probe the time 
dependence of the fluctuations is provided by dynamic 
Rayleigh light scattering. Dynamic light scattering ex- 
periments in fluid layers with negative Rayleigh numbers 
have shown the existence of two modes: a heat mode 
with a decay rate equal to Dxq^ associated with tem- 
perature fluctuations, and a viscous mode with a decay 
rate equal to vq^ associated with transverse momentum 
fluctuations [jjllJI- Thus, for large q where gravity and 
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boundary effects are negligible the coupling between the 
heat mode and the viscous mode causes an enhancement 
of the amplitude of nonequilibrium fluctuations, but does 
not affect the decay rates, in accordance with the original 
predictions of Kirkpatrick ei aZ. 1]. However, for smaller 
q (corresponding to macroscopic wavelengths), gravity 
and boundary effects induce a coupling between the de- 
cay rates of the viscous and heat modes, as originally sug- 
gested by Lekkerkerker and Boon For certain neg- 
ative values of the Rayleigh number the modes can even 
become propagating On the other hand, 

for positive R near the RB instability the nonequilibrium 
structure factor is dominated by a very slow mode with a 
decay rate which vanishes as i? — > i?c 0, 0, E| ■ 

The phenomenon of critical slowing down of the fluctu- 
ations as the RB instability is approached from below has 
been observed experimentally, but the results obtained 
so far are qualitative. Using forced Rayleigh scattering 
AUain et al. poj measured the decay of imposed horizon- 
tal spatially periodic temperature profiles with various 
wave numbers q. The decay time of these imposed de- 
viations from the steady state increased as the temper- 
ature gradient approached the value associated with Re- 
Furthermore, the decay times became larger for periodic 
temperature profiles with q closer to qc- Critical slowing 
down of nonequilibrium fluctuations has also been ob- 
served by Sawada with an acoustic method [23|. In this 
experiment, first a deterministic RB pattern was estab- 
lished and then the cell was turned over so as to change 
the sign of R. The experimental observations were inter- 
preted with an amplitude equation that did not include 
any wave-number dependence. The single decay time ex- 
tracted from the data did indeed increase as the onset 
of convection was approached, but numerical agreement 
with the amplitude equation was poor. Using neutron 
scattering, Riste and coworkers ^2t>2M observed the crit- 
ical slowing down of nonequilibrium fluctuations close to 
the onset of convection in a liquid crystal. Quantitative 
experimental studies have been performed showing the 
slowing down of the dynamics of deterministic patterns 
as the RB instability is approached from above [24L l25l| . 

In the present paper we present a quantitative the- 
oretical and experimental study of the decay rates of 
fluctuations for positive R but below the RB instabil- 
ity. The results of both experiment and theory reveal 
the critical slowing down as the RB instability is ap- 
proached from below. Theoretically, the decay rates can 
be determined, in principle, from a linear stability analy- 
sis of the deterministic Oberbeck-Boussinesq (OB) equa- 
tion s ag ainst perturbations, as has been done tradition- 
ally [l5[l2ll27Ll2ll2ll30t. However, here we derive these 
decay rates by solving the linearized stochastic OB equa- 
tions, obtained by supplementing the deterministic OB 
equations with random dissipative fluxes 0, 0, ^| . 
Because of Onsager's regression hypothesis, the two pro- 
cedures should yield the same result, and indeed they do. 
However, it is important to note that a solution of the 
stochastic OB equations is required to obtain the cor- 



rect amplitudes of the nonequilibrium fluctuations. For 
instance, it turns out that the wave number gmax corre- 
sponding to the maximum enhancement of the intensity 
of the fluctuations as a function of q differs, at least for 
R < Rc, from the value of q at which the decay rate has 
a minimum U] . 

We present in this paper the results of experimental 
shadowgraph measurements 33] of the fluctuations in a 
thin horizontal layer of sulfur hexafluoride heated from 
below but with R < Rc. We show that the decay rates 
of the fluctuations can be obtained from an analysis of 
images acquired at various exposure times, and compare 
the experimental results with the theoretical predictions. 
The data as well as the theory reveal the critical slow- 
ing down of the fluctuations as the onset of convection 
is approached. To our knowledge, previous shadowgraph 
experiments were used primarily to obtain the fluctuation 
intensities 0, . Exceptions are studies of the dynam- 
ics of deterministic patterns above the RB instability that 
used image sequences at constant time intervals [33, |3^ . 
Measurements of the structure factor with the shadow- 
graph method can probe fluctuations at the small wave 
numbers which reveal the influence of the finite geome- 
try. These wave numbers generally are not accessible to 
light scattering. 

We shall proceed as follows. In Sect. ^1 we derive the 
decay rates of the nonequilibrium fluctuations in a fluid 
layer between two horizontal rigid boundaries. To obtain 
analytic expressions we determine the dynamic structure 
factor in a flrst-order Galerkin approximation. This sec- 
tion is a sequel to a previous publication in which 
the static nonequilibrium structure factor was consid- 
ered. In Sect, mil we then use the dynamic structure 
factor to calculate the dependence of shadowgraph sig- 
nals on the exposure-time interval r. In Sect. IIVI we 
describe the experimental conditions and procedures, in- 
cluding the collection and analysis of the shadowgraph 
images. Specifically, we show how the decay rate ^-{q) 
of the fluctuations below the RB instability can be de- 
termined experimentally from shadowgraph images ob- 
tained with two different exposure times. A comparison 
of the experimental results obtained for T-{q) with our 
theoretical predictions is presented in Sect. Finally, 
our conclusions are summarized in Sect. IVTl 



II. DYNAMIC STRUCTURE FACTOR AND 
DECAY RATES OF FLUCTUATIONS 

As was done previously 0, H Hi determine the 
nonequilibrium structure factor by solving the linearized 
stochastic OB equations for the temperature and veloc- 
ity fluctuations. It should be noted that the use of the 
linearized OB equations implies that the critical slow- 
ing down of the fluctuations near the RB instability is 
treated in a mean-field approximation in which the de- 
cay rate will vanish at i? = i?c when the wave number q 
equals a critical wave number gc ■ It is known from the- 
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ory [ToIIbtI and experiment |3| that nonlinear effects will 
cause a saturation of the fluctuating amplitudes when 
they become very large. The effect of nonlinear terms on 
the decay rates depends on the nature of the bifurcation 
in the presence of fluctuations. In the RB case the fluc- 
tuations change the bifurcation, which is supercritical in 
the absence of fluctuations 38] , to subcritical 0, ^3 ■ 
In that case one would expect the decay rates to remain 
finite at the bifurcation; however, this issue is beyond the 
scope of the present paper. 

For the case of a fluid layer between two plates with 
stress- free (i.e., slip) boundary conditions the intensity 
and the decay rates of the nonequilibrium fluctuations 
have been obtained in a previous publication 0. The 
advantage of stress-free boundary conditions is that they 
permit an exact analytic solution of the problem. Here, 
instead, we consider the more realistic case of a fluid 
layer between two rigid boundaries, corresponding to 
two perfectly conducting walls for the temperature and 
with no-slip boundary conditions for the local fluid ve- 
locity. While this case does not permit an exact an- 
alytic solution, a good approximation can be obtained 
by representing the dependence of the temperature and 
velocity fluctuations upon the vertical z coordinate by 
Galerkin polynomials PO]. Specifically, two of us have 
determined the static nonequilibrium structure factor in 
a first-order Galerkin approximation 4]. A first-order 
Galerkin-polynomial solution appears to provide a good 
approximation below the convection threshold for the 
total intensity of the fluctuations 39J. We thus ex- 
tend here the first-order Galerkin treatment considered 
in Ref. in the expectation that it will also provide a 
good approximation for the decay rates and the actual 
amplitudes of the two coupled hydrodynamic modes. 

The decay rates of the hydrodynamic modes can be 
readily obtained by finding the roots in u of the determi- 
nant of the matrix 7i(ijJ, g), defined by Eq. (15) in Ref. 0. 
There are two coupled hydrodynamic modes whose decay 
rates V±{q) are given by 



r±(g)-^(9 +i«)|tftwtw+ 



(2) 



±1 



o-(g4 + 24g2 + 504) 
(g2 + 10)(g2 + i2) 



- 1 



7(g2 + i2)(g- 



f-f2. 



10)2 j 



where cr = v jn, is the Prandtl number, and q — qd is 
the dimensionless horizontal wave number of the fluctu- 
ations. We note that in the previous publication the 
symbol g|| was used to emphasize that it represents the 
magnitude of a two-dimensional wave vector q in the hor- 
izontal Xy-plane. Here we drop the subscript ||, since in 
the present paper the wave vectors will always be two- 



dimensional in the horizontal plane. In Eq. ^-{o) 
represents the decay rate of a slower heat-like mode which 
approaches DtQ^ for large values of q, while r-|_(g) repre- 
sents the decay rate of a faster viscous mode approaching 
lyq^ for large q. The advantage of the Galerkin approx- 
imation is that one can specify the decay rates T±(q) 
explicitly as a function of q. We note that in the first- 
order Galerkin approximation considered here, we find 
only two decay rates, instead of a series of decay rates 
that would be obtained if higher orders were considered 
in the Galerkin expansion. For studying the situation be- 
low the RB instability, where fluctuations decay to zero, 
consideration of the two primary decay rates will be ade- 
quate. As discussed more in detail later, for these Eq. Q 
is a good approximation. 

We are interested in the time correlation function for 
the density fluctuations at constant pressure which is 
directly related to the time-dependent autocorrelation 
function {5T*{c^^,t) ■ 5T{(i' , z' ,t')) of the temperature 
fluctuations by 



CL2 
1 



{5T*{q,z,t)-5T{(^,z',t')) = 

F{q,z,z',At)^^S{ci-c^). 

Here At = \t — t'\ and p is the average fluid density. The 
time correlation function F{q, z, z'. At) is just the inverse 
frequency Fourier transform of the dynamic structure fac- 
tor S{uj, q, z, z') defined by Ecl (9) in Ref. Following 
the steps described in Ref. Q we find after some long 
algebraic calculations that the time correlation function 
-F(<7, z, z'. At) can be expressed as the sum of two expo- 
nentials: 



F(g,z,z',Ai) 
+ i_(g)cxp[-r_(g) A<]} 



^{i+(g)exphr+(g) At] (4) 



z z 



z 



Here the coefficient represents the intensity of the 
fluctuations of the fluid in a local thermodynamic equi- 
librium state at the average temperature T = T (see 
Eq. (2) in Ref. Q). The Galerkin polynomials appear 
in Eq. Q because they have been used to represent the 
dependence of the temperature fluctuations on the verti- 
cal coordinate z. We introduce the dimensionless decay 
rates 



f±(g) = tv r±(g) 



(5) 



where = d"^ / Dt is the vertical thermal relaxation time. 
The dimensionless amplitudes A± (q) in Eq. are then 
given by 
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(g2 + 10) 



i±(g) = ±30- 



24g2 + 504 
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27{q'^ + 24g2 + 594)^2 
28(<72 + 12)2 



r±(g) [r2^(g)-r2_(g)] 



(6) 



In Eq. © , 5*^2 denotes the strength of the enhancement 
of the static nonequihbrium structure factor of the fluid 
in the absence of any boundary conditions: 

^SE=-i^+^^^(VTo)^ (7) 

Here cp is the isobaric specific heat per unit mass H. In 
Eq. 0, as everywhere else in the present paper, aU ther- 
mophysical properties are to be evaluated at the average 
temperature T = T. We note that the amplitudes A±{q) 
depend on the Prandtl number a and on the Rayleigh 
number R not only explicitly in accordance with Eq. lO, 
but also implicitly through expression (0) for the decay 
rates r±(g). 

As noted in the Introduction, the decay rates of the 
fluctuations can also be obtained from an analysis of de- 
viations from steady state on the basis of the determinis- 
tic OB equations, i.e., the OB equations without random 
noise terms. Hence, Eq. ^ for r±(g) is implicit in the 
standard calculation of the convection threshold within 
the same Galerkin approximation employed here |3Clj| . 
However, to obtain the correct amplitudes A±{q), it is 
necessary to solve the stochastic OB equations for the 
fluctuating fields. 

Expression (01 for the decay rates, as well as Eq. © for 
the amplitudes, are rather complicated. They simplify 
considerably for a —> 00. In that limit the decay rate of 
the slower mode r_(g) reduces to 

while the decay rate r+(g) of the faster mode becomes 
proportional to a and is so large that the first exponential 
term in Eq. Q can be neglected. The amplitude A- (q) 
of the remaining exponential term reduces, in the same 
cr — > 00 limit, to a simpler expression to be used later, in 
Eq. (|16() . It is interesting to note that the limit for large 
Prandtl numbers is approached rather fast. For instance, 
at cr = 15, the difference between the actual r_(g) given 
by Eq. Q and the value deduced from the asymptotic 
expression (jSJ is always smaller than 3% for any value of 
q, with the larger deviations of about 3% at q values close 
to (jc. From Eq. ((HJ we find that r_(g) approaches 10 
when q 0, independently of the values of the Rayleigh 
or Prandtl number. The value 10 obtained from our first- 
order Galerkin approximation has to be compared with 
~ 9.87, obtained from the exact theory for the limiting 
value at g — > of this decay rate [33 ■ 



The analysis in the previous paper jj| was specifi- 
cally devoted to the static structure factor, Siq,z,z'), 
which may be obtained by setting At = in Eq. Q 
or, equivalently, by integrating the dynamic structure 
factor S{uj,q, z') over the entire range of frequencies 
u! 4]. However, as discussed in more detail in Ref. 4j, the 
structure factor S{q) that is actually measured in small- 
angle light-scattering or in zero-collecting-time shadow- 
graph experiments is the one obtained after integration 
of the full static structure factor over z and z' [13, so 
that: 

S{q) = -, [ dz [ dz'F{q,z,z\0) (9) 
a Jo Jo 

= ll {A+iq) + i_(g)} = 1^ + 5°EAo^(g)| , 

where A^(q') represents the normalized enhancement of 
the intensity of nonequihbrium fluctuations in the first- 
order Galerkin approximation, as given by Eq. (20) in 
Ref. 0. We note that, in accordance with Eq. @, the 
static structure factor S{q) is expressed as the sum of an 
equilibrium and a nonequihbrium contribution. However, 
the dynamic structure factor and its equivalent, the time 
dependent correlation function given by Eq. can no 
longer be written as a sum of equilibrium and nonequi- 
hbrium contributions. This difference between the static 
and the dynamic structure factor results from the cou- 
pling between hydrodynamic modes due to the presence 
of gravity and boundaries. For the same reason, the 
nonequihbrium intensity enhancement no longer ap- 
pears as a simple multiplicative factor in the expression 
for the dynamic structure factor. The observation that 
the nonequihbrium dynamic structure factor is no longer 
the sum of an equilibrium and a nonequihbrium contri- 
bution already pertains to the dynamic structure factor 
of the "bulk" fluid (i.e., without considering boundary 
conditions), where mixing of the modes is still caused by 
gravity effects 



III. APPLICATION TO THE DEPENDENCE OF 
SHADOWGRAPH SIGNALS ON EXPOSURE 
TIME 

The shadowgraph method provides a powerful tool 
for visualizing flow patterns in Rayleigh-Benard convec- 
tion ||3^ |36|, |42| . Of special interest for the present pa- 
per is that the shadowgraph method can also be used to 
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measure fluctuations in quiescent nonequilibrium fluids 
at very smaU wave numbers [a, [l^, yj, Hfl, Hj] . 

In shadowgraph experiments an extended uniform 
monochromatic Hght source is empfoyed to illuminate the 
fluid layer. Many shadowgraph images of a plane perpen- 
dicular to the temperature gradient are obtained with a 
charge-coupled-device (CCD) detector that registers the 
spatial intensity distribution /(x, r) as a function of the 
two-dimensional position in the imaging plane x and the 
nonzero exposure time r used by the detector to average 
photons for a single picture. An effective shadowgraph 
signal X(x, t) is then deflned as: 



I(x,r) 



/(x, r) - /o(x,t) 
r) 



(10) 



where /o(x, t) is a blank intensity distribution in the ab- 
sence of any thermally excited fluctuations. In practice, 
lo (x, r) is calculated as an average over many original 
shadowgraph images, so that fluctuation effects cancel. 

From a series of experimental shadowgraph signals 
2r(x, t) , the experimental shadowgraph structure factor 
S's(q, t) is deflned as the modulus squared of the two- 
dimensional Fourier transform of the shadowgraph sig- 
nals, averaged over all the signals in the series: S's(q, t) = 
(|Z(q, r)|^). The physical meaning of the shadowgraph 
structure factor iS's(q, r) in the past has been based on 
the assumption that the shadowgraph images are taken 
instantaneously. In the limit t —^ 5*8 (q,r) has been 
related to the static structure factor of the fluid S'(q) by 
a relation of the form i, El S IS 113 



S,{q,0) = T{q) S{q). 



(11) 



In Eq. (|ll|l . T{q) is an optical transfer function that con- 
tains various properties, such as the wave number of the 
incident light, the temperature derivative of the refrac- 
tive index of the fluid, and details of the experimental 
optical arrangement. For the present work a speciflca- 
tion of T(g) is not needed since it will be eliminated in 
the treatment of the experimental data in Section IIVI 
The quantity S (q) in Eq. 1)11(1 equals the static structure 
factor of the nonequilibrium fluid, as defined in Eq. @. 

To determine the dependence of the experimental 
shadowgraph structure factor Ss{q,T) on the exposure 
time r we need to extend the physical optics treatment 
performed by Trainoff and Cannell (4^]. We then find 
that Eq. 1)11(1 is to be generalized to: 



5s(q,T) =r(g) S{q,T), 



(12) 



where T(g) is the same optical transfer function as in 
Eq. (|ll(l and where S{q, r) is a new exposure-time- 
dependent structure factor, which is related to the au- 
tocorrelation function of temperature fluctuations by 

^^5(g,r)5(q-q') = ;^ (13) 

pT pT nd pd 

X dt dt' dz dz' {ST*{q,z,t) ■ ST{q',z',t')). 
Jo Jo Jo Jo 




FIG. 1: Difference SS{q,T)/ SeS^e (see Eg. \Wl below) be- 
tween the theoretical structure factors at i? = 1371 (Eq. 1141 '). 
and at i? = (Eq. lIT^ 'l as a function of q, for three different 
collecting times. The solid curve corresponds to r = ms, the 
dashed curve to r = 200 ms, and the dotted curve to r = 500 
ms. The Prandtl number is u = 34 and tv = 0.74 s. 



In principle, S{q,T) can be evaluated by substituting 
Eq. ||2J) with Eq. Q for the correlation function of the 
density fluctuations into the right-hand side of Eq. I|13() . 
However, for the experimental results to be presented, 
it is sufficient to evaluate S{q,T) for large values of the 
Prandtl number a, in which case the exponential contri- 
bution with decay rate r+ (q) in Eq. I^J can be neglected, 
as discussed before. Retaining only the contribution from 
the slower mode with amplitude A^{q) in Eq. Q and 
performing the integrations, we deduce from Eq. ((13(1 



S{q,f) ^ S^A^iq) 



f_(g)-l + exp[-f f_(g)] 



18f2 rliq) 



(14) 



where f is a dimensionless exposure time, defined by 
f = r/tv with iv the vertical relaxation time introduced 
in Eq. In the hmit t ^ 0, Eq. ifT^ reduces to 

S{q) = S'EA-(g)/36, so that it equals the static struc- 
ture factor measured in small-angle light-scattering or in 
shadowgraph experiments in the zero collecting time ap- 
proximation, as given by Eq. (Q. It is interesting to note 
that, due to the nonzero collecting time, r, even in ther- 
mal equilibrium (i? = 0) the shadowgraph measurements 
present some "structure". From Eq. ((14(1 . we find that 
this equilibrium structure, as a function of the dimen- 
sionless collecting time, f, is given by 



SE{q,f) = S'e — 



5 f{q^ + 10) - 1 + exp[-f (g2 + 10)] 



3f2 



(g2 + 10)2 



(15) 

It can be readily checked that, in the limit r 0, 
Eq. 1(15(1 reduces to the structureless constant (5/6)5e, in 
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agreement with Eq. This value is 17% lower than the 
actual value 5'e, due to the use of a first-order Galerkin 
approximation 

To gain insight into the effect of the exposure time on 
shadowgraph measurements, we show in Fig. ^ the dif- 
ference between the nonequilibrium structure factor as 
obtained from Eq. (|14|) and the equilibrium structure fac- 
tor {R = 0) as obtained from Eq. (|15|l . as a function of 
q, for three different collecting times. Corresponding to 
some of the experimental results to be discussed below, 
we used R — 1371, a vertical relaxation time — (P /Dx 
of 0.74 s, and a Prandtl number a equal to 34. We eval- 
uated the difference for t = 0, t — 200, and r = 500 
ms. The structure factors S{q, r) and S^iq, t) have been 
normalized by dividing them by the product S'eS'^e- 

From Fig.^ we arrive at the following conclusions: (i) 
The main effect of a nonzero collecting time is to lower 
the height of the measured S{q, r); this is expected since 
fluctuations cancel out when larger exposure times are 
used, (ii) An additional effect of a nonzero collecting 
time, which can be observed in Fig. ^ is a displacement 
of the maximum in S{q^ r) to lower q values. This effect 
is mainly due to the subtraction of the equilibrium struc- 
ture given by Eq. H15|l . which decreases with increasing 
q. (iii) Another interesting feature, we infer from Fig. ^ 
is that the dependence on for low q is preserved, while 
the dependence on q~'^ for large q is destroyed for the 
effective structure factor. Actually, at large q, there ex- 
ists a crossover from a g"* dependence to a q~^ depen- 
dence, for nonzero values of r, as will be shown below in 
Eq. ini). 

A simpler formula for S'(g, f ) can be obtained by in- 
troducing into Eq. (|14|l the limiting value of A-i^q) for 
large Prandtl numbers. We then obtain 

^S{q,T) ^ S{q,f) - SE{q,f) 
Se Sj!fE 

5 ft-{q) - l + exp[-ff_(g)] 



3a 



r2 f^(^) 
27g2 



28(g2 + 10) (g4 + 24g2 + 504) - 27 Rq^ 



(16) 



where for f-((?) Eq. ||SJ) for large tr should be used. As 
discussed after Eq. ©, for Prandtl numbers near 15 
the asymptotic limit IjlGI) is already closely approached. 
Equation 1)16(1 reduces, in the limit r — > 0, to the expres- 
sions obtained in Ref. 0. It is interesting to note that, 
in the limit q — > oo, Eq. (|16|l reduces to: 



SS{q,T) q- 
SeSI 



45 1 



NE 



28fCT St" 



(17) 



showing, as mentioned above, that a nonzero collecting 
time changes the asymptotic behavior from to q~^ . 
Note that Eq. H17|l is evidently not valid when r = in 
which case one needs the asymptotic expression given by 
Eq. (30) in Ref. 4]. 



IV. EXPERIMENTS 



A. Apparatus 



To visualize the thermal fluctuations with the shadow- 
graph method one needs to perform the experiments with 
a fluid in which the thermal noise will be large 0, ^3 ■ 
This goal can be accomplished by selecting a fluid in the 
vicinity of its critical point [l^ ITS l33l | . The measure- 
ments reported here were obtained for sulfur hexafluo- 
ride. The apparatus and experimental procedures have 
been described in detail elsewhere ^ . Here we describe 
only those aspects which are specific to the experiment 
with a fluid near its critical point. 

For the details of the cell construction we refer to Fig. 
9 of Ref. 'ss']. Initially we used a diamond-machined 
aluminum bottom plate which could be positioned with 
piezo-electric elements. The bottom-plate thermistors 
were embedded approximately 0.64 cm below its top sur- 
face. Even though to the naked eye this plate had a near- 
perfect mirror finish, the tool marks from the diamond 
machining imposed a preferred direction on the fluctua- 
tions below the onset of convection. Thus we placed an 
optically flat sapphire of thickness 0.318 cm on top of the 
aluminum plate. A thin silver fllm was evaporated on the 
top surface of this sapphire to provide a mirror for the 
shadowgraphy. 

Initially we used a sapphire of thickness 0.952 cm for 
the cell top. In addition a sapphire of thickness 1.90 cm 
was in the optical path of the shadowgraphy and provided 
the top window of the pressure vessel. In combination 
with the very small cell thickness used in these experi- 
ments, it turned out that the optical anisotropy of these 
randomly oriented sapphires introduced an anisotropy of 
the shadowgraph images which obscured the rotational 
symmetry of the fluctuations. In order to minimize this 
effect, we replaced the pressure window by a fused-quartz 
window and used an optically flat sapphire of 0.318 cm 
thickness for the cell top. Under these conditions we 
found that the fluctuations below the onset of convection 
yielded a structure factor which was nearly invariant un- 
der rotation, as can be seen from Fig. |S1 below. 

The cell spacing was fixed by a porous paper sidewall 
with an inner (outer) diameter of 2.5 (3.5) cm. Since the 
top sapphire was supported along its perimeter which 
had a diameter of 10 cm (i.e., considerably larger than 
the cell wall), the force exerted on the cell top by the 
bottom plate and the wall caused a slight bowing of the 
initially flat top. Over the entire sample diameter this 
yielded a radial cell-spacing variation corresponding to 
about one circular fringe when illuminated with an ex- 
panded parallel He-Ne laser beam. This variation of the 
thickness by about 0.3 ^m assured that convection would 
start in the cell center, rather than being nucleated in- 
homogeneously near the cell wall. Assuming a parabolic 
radial profile for the cell spacing, we estimate that the 
spacing was uniform to much better than 0.1% over the 
1.3x1.3 cm^ area near the cell center which was actu- 
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FIG. 2: The temperature-density plane near the critical point 
of SFe. The dashed line indicates the coexistence curve sep- 
arating liquid and vapor. The vertical dotted line is the 
critical isochore. The solid circle is the critical point with 
Tc = 45.567° C, Pc = 37.545 bar, pc = 0.742 g/cm^ The 
solid line represents the isobar P = 38.325 bar used in our 
measurements, and the heavy section of this line shows the 
temperature-density range on the isobar which is spanned by 
the sample with thickness d — 34.3 fim when AT = 0.22° C. 



ally used for the shadowgraph images. The actual sam- 
ple thickness was measured interferometrically [s^ and 
found to be 34.3 ^m. 



B. Properties of SFe near the critical point 

The thermodynamic properties of sulfur hexafluoride 
in the critical region can be calculated from a equation 
developed by Wyczalkowska and Sengers For the 

viscosity we used a fit of data from Refs. |44| and 
to a smooth function. This approach neglects a small 
anomaly of the viscosity at the critical point. We used 
fits of smooth functions to the conductivity data from 

Refs. miiiiiiiiiiSiil. 

The measurements were made at constant pressure P 
and at constant mean sample temperature T. The mean 
temperature was adjusted so that the density p{T) was 
the critical density pc = 742 kg m""^. The imposed tem- 
perature difference AT caused a density variation of the 
sample along an isobar. This is illustrated in Fig. |21 
for the conditions of the present experiments, namely 
for P = 38.325 bar and f = 46.50°C. At p = p^, 
we have a = 34.0, Dt = 1.59 x 10~^ cm^/s, and 
U = (f/Dr = 0.738 s. 




46.2 46.4 46.6 46.8 



T (°C) 

FIG. 3: The thermal expansion coefficient a along the isobar 
P = 38.325 bars. The heavy part of the line terminated by 
two open circles indicates the temperature range spanned by 
the sample with d = 34.3 pm and AT = ATc = 0.44° C. 



C. Symmetric departures from the 
Oberbeck-Boussinesq approximation 

Much of the theoretical work on RBC has been done in 
the OB approximation [2fil f27l | , which assumes that the 
fluid properties do not vary over the imposed tempera- 
ture interval, except for the density where it provides the 
buoyant force j^^. Non-OB effects have been considered 
by a number of investigators, and most systematically 
by Busse [S^. At leading order, they break the reflection 
symmetry of the system about the horizontal midplane, 
and at the onset of convection they yield a transcriti- 
cal bifurcation to a hexagonal pattern 1^3, instead of 
the roll pattern of pure OB convection. When the mean 
temperature corresponds to the critical isochore, p = pc, 
this effect is of modest size. Although several proper- 
tics contribute, we illustrate this by showing the isobaric 
thermal expansion coefhcient in Fig. |31 along the isobar 
of Fig. 121 as an example. One can approximate a as a 
sum of two contributions, one of which is anti-symmetric 
and the other one symmetric about the mean tempera- 
ture (and thus approximately also about the horizontal 
midplane) of the sample. Only the anti-symmetric part 
is considered in the theory [53|. Its smallness near the 
onset of convection is seen from the similar values of a 
at the top and bottom of the cell (open circles in the fig- 
ure) . A quantitative calculation of the parameter V in- 
troduced by Busse IS^ (see Eqs. (13) in Ref. yields 
V = —0.23. This value indicates that, within our exper- 
imental resolution, only rolls should be seen near onset. 
This is indeed the case However, the variation of a 
which is symmetric about the midplane is quite large. It 
does not break the reflection symmetry about that plane 
and thus permits the existence of rolls at onset rather 
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FIG. 4: The Rayleigh number R (see Eq. (0) along the isobar 
P — 38.325 bar for d — 34.3 /im. From top to bottom, the 
curves are for AT/ATc = 1.000, 0.860, 0.699, 0.430, and 0.269 
with ATc = 0.44° C. 

than requiring a hexagonal plan form. Behavior simi- 
lar to that of a is found for the specific heat cp and for 
the Rayleigh number. At present there is no theoretical 
treatment of these higher-order non-OB effects. Thus we 
proceed empirically by examining the variation of R as 
a function of vertical position or local sample tempera- 
ture. If we neglect the temperature dependence of the 
thermal conductivity and assume that the local temper- 
ature in the fluid layer still varies linearly as a function 
of z, we obtain the Rayleigh number profiles shown in 
Fig. 13 The curve at the top is for the experimental 
value AT = ATc = 0.44°C corresponding to the onset 
of convection. The other curves, from top to bottom, 
correspond to AT/ATc =0.860, 0.699, 0.430, and 0.269. 
For AT = 0.44°C one sees that the local R{pc) ~ 3190 
far exceeds the value Rc = 1708 for the uniform sys- 
tem. On the other hand, near the top and bottom of the 
sample the local R is well below the onset of convection 
for the homogeneous system. The data in Fig. |S1 show 
that, though we find rolls above threshold, non-OB ef- 
fects cannot be completely neglected in our experiments. 
We shall return to the influence of non-OB effects on the 
comparison between experiment and theory in Sect. IVl 



D. Sample temperature 

The bath temperature Tbath and the bottom-plate tem- 
perature Tbp were adjusted so as to hold the mean sam- 
ple temperature constant. This temperature was chosen 
so that the mean density corresponded to the critical 
density pc- Because of the small sample thickness the 
thermal resistance of the sample was comparable to that 
of the top and bottom confining plates. This fact re- 




FIG. 5: Shadowgraph signals (left column) and the mod- 
uli squared of their Fourier transforms (right colum) for 
AT = 0.189 (top row) and AT = 0.378 K (bottom row). 
The exposure time was 500 ms. 



quired a special procedure to assure that the sample was 
indeed at the temperature corresponding to pc- Before a 
run at a given fixed pressure was started, we measured 
the power of shadowgraph images of the fluctuations at 
a fixed imposed AToxt = (Tbp — Tbath) as a function of 
the mean temperature Text = (Tbp + Tbath)/2 of the sys- 
tem consisting of the bottom plates, the sample, and the 
top plate as determined by the bath temperature and the 
bottom-plate thermistor. This temperature difference is 
the sum of those across the bottom aluminum plate ATai, 
across the boundary between the aluminum plate and the 
bottom sapphire ATbi, across the bottom sapphire ATgb, 
across the sample AT, across the top sapphire ATgt, and 
across a boundary layer above the top sapphire in the wa- 
ter bath ATb2 . From estimates of the thermal resistances 
of these sections we find AT/ AToxt — 0.473. Thus, at the 
onset of convection, we measured AToxt. c = 0.930°C and 
deduced ATo = 0.44°C. 



E. Analysis of shadowgraph images 

At each AT, three shadowgraph-image sequences 
Iz{x,Tj), i = 1,...,N with N = 1024, and j = 0,1,2 
were acquired. For each sequence, a different exposure 
time Tj was used, namely tq = 0.500 s, ti = 0.350 s, 
and T2 = 0.200 s. The time interval St between the im- 
ages was typically one or two seconds, which was large 
enough for the images to be nearly uncorrelated. The 
images of each sequence were averaged to provide a back- 
ground image /o(x, Tj), as discussed in Sect. IIIII Then, 
for each image of the sequence, a dimensionless shadow- 



9 



graph signal Ti(x, tj) was computed, in accordance with 
Eq. IjlUI) . The mean (over x) value of a typical 2^ (x, tj ) 
was within the range ±0.01, indicating adequate stability 
of the light intensity and image-acquisition system. Next, 
the two-dimensional Fourier transform of each shadow- 
graph signal in the sequence was computed, and the mod- 
ulus square calculated, obtaining series of |Zi(q, rj)p for 
further analysis. Typical shadowgraph signals and the 
squares of the moduli of their Fourier transforms are 
shown in Fig.jS] The nearly- uniform angular distribution 
of the transforms illustrates the rotational invariance of 
the Rayleigh-Benard system. 

As explained in Sect. IIIII due to the rotational in- 
variance of RB convection in the horizontal plane, the 
modulus squared Fourier transformed shadowgraph sig- 
nals have rotational symmetry, and for an infinitely ex- 
tended sample they would depend only on the modulus 
q of the wave vector q. However, the finite spatial extent 
of the images leads to random angular fluctuations of 
|Xj(q, rj)p. To reduce these fluctuations, we performed 
azimuthal averages |Xi(q, Tj)p over thin rings in Fourier 
space (the angular average is denoted by the overline and 
depends only on q and Tj). Now for each shadowgraph 
measurement Ii{yi,Tj) the integral 

P,{Tj)^ 27Tq |I,(q,Tj)|2 dq (18) 
Jo 

is the total power and, by Parseval's theorem, has to 
be equal to the variance of the original Zi(x, Tj). We 
used Eq. (|18|l as a check of consistency for the entire 
procedure of taking the Fourier transform, of calculating 
the modulus squared and the azimuthal average over thin 
rings, and of assigning to each ring a q value. Finally, 
we averaged over the N individual |Ii(q, Tj)\'^ of each tj 
series, to obtain the experimental shadowgraph structure 
factor Ss{q,Tj) = (|Ii(q, r,)|2). 

For the remainder of this paper we shall use the exper- 
imentally determined cell spacing d — 34.3/im to scale 
q according to q — qd. Two examples of the product 
2TrqSs{q,Tj), both for P = 38.325 bar, f = 46.5°C, 
and AT = 0.189°C, are shown in Fig. The sohd 
squares are for T2 = 0.200 s, whereas the solid circles are 
for To = 0.500 s. As discussed in the previous section, 
the relationship (|12|l between the experimental shadow- 
graph structure factor Ss{q,Tj) and the corresponding 
fluid structure factor S{q, Tj) involves the optical transfer 
function T{q). Nevertheless, we observe some qualitative 
agreement of the experimental results for 2'KqSs{q,T) in 
Fig. IHl and the theoretical results for S{q, r) shown in 
Fig. ^ The expected maximum near q = qc \s present, 
and the decrease as q vanishes is consistent with the pre- 
dicted (f' dependence. The longer averaging of the ran- 
dom fluctuations diminishes 27rg Ss{q, t), which was one 
of the features noted after Fig. Q for S{q, r). 

The experimental structure factor Ss{q,Tj, AT) de- 
pends on the Rayleigh number and, hence, on AT. For 
AT = 0, instrumental (mostly camera) noise is expected 
to be the dominant contribution, i.e., the contribution 
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FIG. 6: Experimental shadowgraph structure factor 
2TTqSs{q,T) as a function of the dimensionless wave num- 
ber q. These results are for P = 38.325 bar, f = 46.5° C, 
and AT = 0.189°C. The solid squares (solid circles) are for 
r — 0.200 s (r = 0.500 s). The open symbols are the corre- 
sponding background measurements for AT — 0. 



from the equilibrium fluctuations is negligible and the 
theoretical equilibrium structure, given by Eq. (|15|l . is 
expected to be unobservable in our experiments. Exam- 
ples of the shadowgraph structure factors for AT = 
are shown by the open symbols in Fig. El One sees that 
2TTqSs{q, Tj, 0) is well represented by a straight line, cor- 
responding to white noise. At large q, 2'!TqSs{q,Tj,0) 
merges smoothly into the data for the 2TrqSs (q, Tj , AT) 
with the same Tj, as one would expect. 

A detailed comparison between the experimental and 
the theoretical structure factors requires knowledge of the 
optical transfer function, T{q), which depends on the de- 
tails of the experimental optical arrangement. It involves, 
for instance, the size of the pinhole and the focal length of 
the lens used to make the "parallel" beam, and the spec- 
tral width of the light source. |4(i] In the present work 
the spatial structures to be determined (the fluctuation 
wavelengths) had length scales [typically 0(50 /Ltm)] that 
are one or two orders of magnitude smaller than those 
of more conventional RB experiments. For this reason, 
we found it difhcult to obtain T(q) for our instrument 
with sufficient accuracy to avoid significant distortion 
of S{q). We circumvented the difficulty of the optical 
transfer function by deriving the dynamic properties of 
the fluctuations from ratios of Ss{q, Tj) with different val- 
ues of Tj. To account for the instrumental white noise, 
before taking such ratios, we subtracted the measured 
Ss {q, Tj , 0) in the absence of a temperature gradient from 
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FIG. 7: The ratio 7^(g, , tq, AT) (see Eq. (|1SJ) based on 
Tj = 0.200 s and tq — 0.500 s. From top to bottom the data 
are for AT = 0.118, 0.189, 0.307 and 0.378°C. The onset of 
convection occurred at ATc = 0.440° C. 



Ss{q,Tj,AT) to yield 

SSsiq, rj,AT) = S,iq, r,, AT) - rj,0). (19) 

After this background subtraction, we formed for j = 1,2 
the ratio 



n{q,T,,To,AT) 



— 3 



t] SSsiq,T,,AT) 
t2 dSMro,^T) 



(20) 



where tq is the longest exposure time tq = 0.500 s. The 
ratios TZ obtained for tj — 0.200 s as a function of q 
are shown in Fig. for four values of the temperature 
difference AT. The shadowgraph transfer function T{q) 
cancels and is no longer contained in TZ. In addition, 
for the ratios TZ, it is irrelevant whether we use the def- 
inition of the structure factor considered in Sect. ^ or 
the shadowgraph definition displayed in Fig. El the lat- 
ter including a factor 27rq. Thus, we are allowed to use 
Eq. lO for SSs, so that 



TZ{q,Tj,To,AT) 



'Tj'^iq) - 1 + L^xp(-Tjr(q)) 



(21) 



ToT{q) - 1 + cxp(-ror(g)) 

One sees that TZ depends only on tq, Tj, q, and T{q, AT). 
At each q and AT, only the decay rate T is unknown and 
thus can be determined from the experimental value of 
TZ. 



V. EXPERIMENTAL RESULTS AND 
COMPARISON WITH THEORY 

In Fig. |S1 we show the decay rate F as a function of 
q. The symbols represent the experimental values de- 



FIG. 8: Results for T{q, AT), in units of inverse seconds, ob- 
tained from the values of TZ shown in Fig. [7| (see Eq. 1211 '). 
The symbols are the same as in Fig. |7| The top (dash-double- 
dotted) curve is the prediction for AT = 0. The remaining 
curves are the theoretical predictions (see Eq. ©) for the 
values of AT of the data (see Fig.0. 



duced from the data for TZ displayed in Fig.[7|by solving 
Eq. H21() . From top to bottom, the data sets are for 
AT/ ATc = 0.269, 0.430, 0.698, and 0.861. The curves 
represent the theoretical values calculated from Eq. © 
with tv = 0.738 s (the value derived from the fluid proper- 
ties at the mean temperature and at the critical density). 
The topmost (dash-double-dotted) line is for equilibrium: 
AT = {i.e., R = 0). The remaining four curves are for 
the values of AT/ ATc of the data sets, if we adopt the 
Boussinesq estimate 



Rob = (AT/ATe)i?c 



(22) 



with i?c = 1730 (the value obtained from the Galerkin 
approximation used in this paper [s^l) for the Rayleigh 
number in Eq. One sees that the predictions, based 
on the Boussinesq approximation, do not agree very well 
with the experiment. However, both theory and exper- 
iment reveal clearly a minimum of F near g ~ 3 which 
becomes more pronounced as AT approaches ATc. For 
the larger values of AT/ATc the experimental values of 
r have a maximum near q — 5. We expect that this is 
due to nonlinear effects in the physical system which lead 
to second-harmonic generation. This phenomenon is not 
contained in the theory. 

We believe that the major disagreement between the 
measurements and the calculation is due to symmetric 
non-OB effects discussed in Sect. IIV CI Since there is 
no quantitative theory, we proceeded empirically and ex- 



11 




2 4 6 

FIG. 9; Results for Tt^ as a function of fqq, using the val- 
ues tv = 0.551 s and fq = 0.944 from the least-squares fit 
described in the text. The data and symbols correspond to 
those in Fig. |H1 The curves indicate the corresponding theo- 
retical results obtained from Eq. (|HJ by using the scale factors 
of the Rayleigh number fa^k from the least-squares fit. 
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FIG. 10: Values of the fit parameter fa^k obtained from the 
least-squares fit described in the text. Solid symbols corre- 
spond to r based on exposure times tq = 0.500 s and tj = 
0.200 s. For the open symbols Tj was 0.350 s. The solid curve 
is the polynomial fa = 2.74 - 2.87Ar/Arc 1.13(Ar/Arc)^ 
which passes through the point (1,1). The dashed curve rep- 
resents {R)/{R)c and was calculated from values of R like 
those shown in Fig. |5] 



plored the possibility that non-OB effects can be accom- 
modated to a large extent by multiplying the experi- 
mental AT/ATc used to estimate R by an adjustable 
scale factor. We introduced an adjustable parameter fji k 
which was allowed to be different for each ATk and used 

^fe = fR,kRoB (23) 

for the Rayleigh number. In addition, we introduced a 
single adjustable scale factor for all data sets which ad- 
justed the length scale of the experiment so as to yield a 
corrected wave number 

Qcoir = fqq- (24) 

to be used in the fit of the theory to the data. We expect 
fq to compensate for experimental errors in the cell spac- 
ing and in the spacing between the pixels of the images, 
to be within a few percent of unity, and to be the same for 
the runs at all AT. Finally, we treated the vertical relax- 
ation time tv as an adjustable parameter. We note that 
tv = d? / Dt depends on the cell spacing, and thus any 
error in the length scale will lead to an error in the time 
scale tv One also might expect tv to depend on AT be- 
cause quadratic non-OB effects would be larger at larger 
AT; but it turned out that a single value for tv for the 
runs at all AT was sufficient to describe all the data. We 



carried out a simultaneous least-squares fit of Eq. H21|l to 
a group of 11 data sets of F, like those in Fig. |S1 based 
on To = 0.500 s and Tj = 0.200 s. Each data set was for 
a different AT^, and collectively they spanned the range 
0.27 < ATfc/ATc < 0.86. A separate fit was done to the 
second group of 9 available data sets based on tq = 0.500 
s and Tj — 0.350 s. For each group we simultaneously 
adjusted fq,t^, and the eleven or nine f^.M- We obtained 
the same result fq = 0.944 ± 0.010 from both groups. 
The fit also gave U = 0.551 ± 0.02 (0.565 ± 0.029) s for 
the group based on Tj — 0.200 (0.350) s. Qualitatively 
consistent with the expected influence of the symmetric 
non-OB effects, the fitted value of tv is somewhat smaller 
than the value 0.738 s estimated for p ^ pc and the ex- 
perimental value d = 34.3 pm. As said above, part of 
this difference is attributable to the error in d indicated 
by the result obtained for fq. In Fig. |51 we show the 
results for the product F = rt„ for the four examples 
displayed in Fig.|Slas a function of fqij, together with the 
corresponding predictions generated by using the values 
of fn^k from the fit. Except for the second-harmonic con- 
tribution at large q and AT, the adjusted theory agrees 
quite well with the data. 

The values obtained for fn^k are given in Fig.^| The 
two groups {Tj = 0.200 s and tj — 0.350 s) agree very 
well with each other, showing that consistent results are 
obtained with different exposure times. Also shown is a 
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fit to the data sets with Tj = 0.200 s in which the indi- 
vidual //j_fc were replaced by a quadratic function which 
was forced to pass through fn^k = 1 at AT — ATc. This 
fit gave 



fR{AT/AT,) = 2.74 - 2.87 



AT 

at; 



1.13 



/AT 

Iai; 



(25) 



This two-parameter representation of fjj fits the data 
equally well. 

One sees that fn is largest at the smallest AT, indicat- 
ing that the estimate R = Rqb (see Eq. (|22ll ') becomes 
worse as AT decreases. This is to be expected because 
the approximation for fji must approach unity as AT 
approaches ATc in order for r(gc, ATc) to vanish. As a 
simplest empirical attempt to include non-OB effects in 
the comparison between experiment and theory, we de- 
fine, in analogy to Eq. H22I) . a non-Boussincsq Rayleigh 
number 



(26) 



where as before Rc = 1730 and where the angular bracket 
indicates an average over the spanned temperature (and 
thus density) range along the isobar. The averaged 
critical Rayleigh number (i?(p))c is equal to {R{p)) for 
AT = ATc. It turns out that {R{p))c = 1120 for our 
experiment. This approximation corresponds to a re- 
defined 



/i?(AT) 



ATc {R{p)) 



AT (Rip)), 



(27) 



We note that by definition /ij(ATc) = 1 and thus 
R — 1730 as it should be. Equation (|?7jl is plotted in 
Fig. (|10|l as a dashed line. One sees that this simplest 
non-OB model accounts very well for the experimental 
data of Jr. Of course, it would be very helpful to have 
a proper theory (rather than an empirical model) of this 
interesting effect. 



VI. SUMMARY 

In this paper we have reported on experimental and 
theoretical studies of the dynamics of thermal fluctua- 



tions below the onset of Rayleigh-Benard convection in 
a thin horizontal fluid layer bounded by two rigid walls 
and heated from below. Starting from the fluctuating 
linearized Boussinesq equations, we derived theoretical 
expressions for the dynamic structure factor and the de- 
cay rates and amplitudes of the hydrodynamic modes 
that characterize the dynamics of the fluctuations. The 
dynamic structure factor is dominated by a slow mode 
with a decay rate that vanishes as the Rayleigh number 
R becomes equal to its critical value Rc for the onset of 
convection. 

We used the shadowgraph method to determine the 
ratio TZ of shadowgraph structure factors obtained with 
different camera exposure times. From the theoretical 
results for the dynamic structure factor, we derived a 
relationship between this ratio and the decay rates of the 
fluctuations. Using this result and the experimental val- 
ues of TZ, we obtained experimental decay-rate data for 
a wide range of temperature gradients below the onset of 
Rayleigh-Benard convection in sulfur hexafluoride near 
its critical point. Quantitative agreement between the 
experimental decay rates and the theoretical prediction 
could be obtained when allowance was made for some 
experimental uncertainty in the small spacing between 
the plates and an empirical estimate was employed for 
symmetric deviations from the Oberbeck-Boussinesq 
approximation which are expected in a fluid with its 
mean density on the critical isochore. 
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